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ABSTRACT 

We have obtained a Gemini South T-ReCS Q a -band (18.3 /j,m) image and a Spitzer MIPS SED- 
mode observation of HD 181327, an F5/F6V member of the ~12 Myr old /3 Pictoris moving group. 
We resolve the disk in thermal-emission for the first time and find that the northern arm of the disk 
is 1.4x brighter than the southern arm. In addition, we detect a broad peak in the combined Spitzer 
IRS and MIPS spectra at 60 -75 fim that may be produced by emission from crystalline water ice. 
We model the IRS and MIPS data using a size distribution of amorphous olivine and water ice grains 
(dn/da oc a -2 25 with a m i n consistent with the minimum blow out size and a max = 20 /im) located 
at a distance of 86.3 AU from the central star, as observed in previously published scattered-light 
images. Since the photo-desorption lifetime for the icy particles is ^1400 yr, significantly less than 
the estimated ~12 Myr age of the system, we hypothesize that we have detected debris that may be 
steadily replenished by collisions among icy Kuiper belt object-like parent bodies in a newly forming 
planetary system. 

Subject headings: stars: individual (HD 181327) — stars: circumstellar matter — planetary systems: 
formation 



1. INTRODUCTION 

Debris disks are dusty, gas-poor disks around main 
sequence stars. The estimated lifetimes for their con- 
stituent grains are typically shorter than the ages of their 
central stars by orders of magnitude, suggesting that the 
observed dust is debris, replenished from a reservoir such 
as by collisions among parent bodies or sublimation of 
comet-like bodies. To date, Infrared Astronomical Satel- 
lite (IRAS) and Spitzer studies have identified hundreds 
of debris disk candidates via thermal-emission from cir- 
cumstellar dust (e.g. Hillenbrand et al. 2008; Rhee et al. 
2007). Approximately 95% of these candidates possess 
dust with cool color temperatures (T gr < 100 K), similar 
to that expected for dust in the solar system's Kuiper 
Belt (Meyer et al. 2007). Broad-band near- infrared pho- 
tometry and low resolution spectroscopy suggest that 
KBOs possess ices on their surfaces (Schaller & Brown 
2007; Delsanti et al. 2006). The composition of dust in 
debris disks may be similar but has not yet been well- 
quantified (Chen et al. 2006). 

High resolution scattered-light and thermal-emission 
imaging and spectroscopy may reveal disk structure and 
provide insight into grain properties via (1) brightness 
peaks in scattered-light or thermal-emission images that 
indicate local dust density enhancements and/or dust 
temperature asymmetries and (2) spectral features at 
mid- and far- infrared wavelengths. Only a handful of 
debris disks have been imaged both in thermal-emission 
and scattered-light thus far: Pictoris, HD 32297, HR 
4796A, HD 107146, and Fomalhaut. There may be tan- 
talizing evidence to suggest that ices are present in these 
debris disks. Highly porous grains, similar to those used 
to model the infrared spectra of comets, have been used 
to model the (3 Pictoris infrared spectrum and spectral 
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energy distribution (Chen et al. 2007; Li & Greenberg 
1998). Excess thermal-emission in the inner region of 
the HD 32297 disk may indicate the presence of ice at 
temperatures near sublimation (Fitzgerald et al. 2007). 
However, moderate resolution, far-infrared spectra are 
needed to detect unambiguously emission from icy grains 
that is expected to be abundant in these systems. 

HD 181327 is an F5/F6V star with a Hipparcos- 
determined distance of 50.6 pc. It has been identified 
as a member of the (3 Pictoris moving group, with an es- 
timated age of ~12 Myr (Zuckerman et al. 2001). Man- 
nings & Barlow (1998) first identified HD 181327 as a 
possible debris disk based on the detection of IRAS 25 
and 60/im excesses. Hubble Space Telescope (HST) NIC- 
MOS and ACS coronagraphic imaging has resolved an in- 
clined ring of dust in scattered-light at a distance of 86.3 
AU with a width of 36 AU and an inclination, i = 31.7° 
from face-on. The scattered light is ^0.17% of the inci- 
dent stellar light and possesses a brightness asymmetry 
that is well fit using an asymmetric scattering parameter, 
g = 0.30±0.03, at 1.1 [im, suggesting that the grains pro- 
ducing the scattered-light are small, with radii smaller 
than 1 /jm (Schenider et al. 2006) if the grains are uni- 
formly distributed in the disk. Spitzer IRS observations 
measure a steeply rising 5-35 (im spectrum that is well 
characterized by an 80 K single-temperature black body 
and a fractional infrared luminosity, Lj^J — 3.1 x 10 -3 
(Chen et al. 2006), commensurate with that measured 
toward (3 Pictoris. We have carried out a high resolu- 
tion multi-wavelength imaging and spectroscopy study 
to infer the grain properties around HD 181327. 

2. OBSERVATIONS 

2.1. T-ReCS Observations 

We used the T-ReCS imager on Gemini-S to observe 
HD 181327 on the night of 2006 July 21 (GS-2006A-Q- 
43) . These observations were interleaved with exposures 
of a Ind, a photometric and point-spread function (PSF) 
reference. Both targets were imaged in the Q a band 
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(Ao = 18.3 /im) with an ABAB chop-nod sequence (15"- 
throw at P.A. = 30°). The total guided, on-source ex- 
posure times were 4865 s and 456 s for HD 181327 and a 
Ind, respectively. 

The data were reduced using custom Python routines. 
In order to suppress sky and instrumental emission, we 
used standard addition and subtraction tehcniques for 
combining the exposures in each AB nod pair. For 
HD 181327, two of the resulting double-difference images 
with abnormally high variances (more than twice the me- 
dian value) were discarded, reducing the on-source expo- 
sure time by 43 s. Each of the remaining nod-pair im- 
ages was then weighted by its variance, and images with 
common pointings were summed. These resulting images 
were then processed to remove electronic offsets. After 
masking a 3"-radius region around the target, the me- 
dian values in each row and (20-column) amplifier chan- 
nel were used to iteratively solve for and remove offsets. 
Each image was then registered at the one-pixel level 
with Gaussian centroids, and the results were shifted 
and summed to a single image. This image was simi- 
larly processed to remove offsets in each row and column 
(as opposed to each row and 20-column channel in the 
earlier step). 

Detailed examination of the resulting image of 
HD 181327 revealed additional electronic pattern noise, 
seen as diagonal banding. The power spectrum of the 
image contained peaks at (nv x ,Vy), where v~ l = 20pix 
(corresponding to the amplifier channel width), Vy 1 = 
55pix, and n — 1,2, • •-. After masking a 3"-radius re- 
gion around the star, Hermite-symmetric Fourier modes 
at spatial frequencies within 1.5 pix -1 of the n = 1-3 
peaks were fit to the image of HD 181327 and subtracted. 
Additionally, all power at frequencies above the diffrac- 
tion limit v t oi = -D/A m i n (taking A m j n = 17.57/xm for 
Q a ) was removed from the images of both HD 181327 
and a Ind. The results of this process are shown in 
Figure [TJ While the subtraction of Fourier modes at 
the power spectrum peaks reduced the 1-tr noise level by 
1.4% (3.3% in smoothed images), the primary purpose of 
this step was to remove spatial correlations in the noise 
which might masquerade as faint circumstellar emission. 

We used aperture photometry of a Ind for photometric 
calibration. We computed the flux zero point using an 
aperture of 3" 6 radius, and a Q a flux value of 6.032 Jy. 
An examination of the enclosed flux as a function of air- 
mass showed no clear trend, suggesting time- variable tel- 
luric extinction. The mean airmasses of the HD 181327 
and a Ind exposures were 1.18 and 1.22, respectively. We 
elected not to apply a correction for the expected differ- 
ence in mean extinction between the two stars. Rather, 
we note that time- varying extinction is likely to dominate 
the systematic photometric errors. We note that aper- 
ture photometry of a Ind in nod-pair difference images 
had a variation of 11%, which we take as an estimate of 
the systematic uncertainty in the photometric zero point. 

2.2. MIPS SED-mode Observations 

We obtained follow-up Spitzer (Werner et al. 2004) 
MIPS (Rieke et al. 2004) SED-mode observations 
(55.0-90.0 /im; A/AA - 20) of HD 181327 (AOR key: 
16171008) on 15 April 2006, using the 19.6" x 157" 
slit/grating and 10 cycles of 10 s integrations. Each ob- 



serving cycle consists of six pairs of 10 s exposures with 
the slit position alternated between the object and blank 
sky 1' away. To remove the background, each sky expo- 
sure is subtracted from the immediately preceding frame 
containing the object. Exposures of an internal calibra- 
tion source are interspersed within the cycle to track the 
varying response of the 70 micron detector array. The 
raw MIPS data were corrected for distortion, registered, 
mosaicked, and flat fielded using the MIPS instrument 
team's data analysis tool (DAT; Gordon et al. 2005). 
The spectrum was extracted using a 5-pixel extraction 
aperture and calibrated using the procedure and calibra- 
tion files described in Lu et al. (2008). The processing 
steps include applying the dispersion solution, an aper- 
ture correction (to account for slit losses) and a flux cal- 
ibration based on ~20 IR standard stars (Gordon et al. 
2007). 

We use these newly obtained MIPS SED-mode data in 
conjunction with previously published Spitzer /IRS ob- 
servations to model the SED of HD 181327. The obser- 
vations and data reduction for the IRS data are described 
in Chen et al. (2006). We estimate the HD 181327 stel- 
lar photosphere using a Kurucz stellar atmosphere model 
with T ff = 6500 K, log g = 4.0, and solar metallicity, nor- 
malized to the first 10 data points (A < 6 pLm) in the IRS 
SL spectrum. We assume a 5% uncertainty in our stellar 
photosphere model. 

3. THERMAL IMAGING RESULTS 

We plot our g a -band images of HD 181327 and PSF 
reference star a Ind in Figure [T] and their enclosed flux 
curves in Figure^, with 1-cr uncertainties derived from 
the background noise. We subtract the PSF star image 
from our HD 181327 image, assuming that the stellar 
point source has a flux of 53 mJy (estimated from our 
stellar photosphere model), revealing thermal-emission 
from the disk within apertures <3" (Figure 3). The lack 
of excess interior to 0"4 suggests that our photospheric 
flux estimate is accurate. We measure an excess Q a - 
band flux of 61.5±2.7mJy within a 3" circular aperture, 
consistent with our analysis of the Spizter /IRS spectrum. 

HST ACS and NICMOS scattered-light imaging re- 
solve an inclined ring of dust around HD 181327 with (1) 
a surface brightness that falls off rapidly interior to the 
86.3 AU ring and (2) an azimuthal brightness asymme- 
try that is well-modeled using a Henyey-Greenstein phase 
function. These observations suggest that the inner re- 
gion of the disk possesses little dust and assume that the 
distribution of dust in the ring is azimuthally symmet- 
ric. In Fig. [3l we indicate the extent of the scattered- 
light ring with ellipses corresponding to the inner and 
outer e -folding radii of the deprojected scattered-light 
profile (|Schneider et al.ll200l §5.2). 

The northern-arm of the disk appears brighter than 
the southern arm, contrary to the equal brightness ex- 
pected from an azimuthally-symmetric dust ring. We 
measure the flux in annular wedges, centered on the mi- 
nor axis, with outer boundaries defined by the outer e- 
folding radius measured in scattered light. We measure 
a fluxes of 18 ± 1 mJy and 13 ± 1 mJy in a 120° wedges 
in the northern and southern arms, respectively, corre- 
sponding to a flux ratio of 1.4 ± 0.1 between the north- 
ern and southern arms. Therefore, the observed large 
scale azimuthal variation in scattered light may not be 
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purely generated by the dust scattering phase function 
as previously assumed. Instead, the dust density and/or 
temperature may be higher in the northern arm than 
in the southern arm. Faint diffuse emission is observed 
in scattered-light with ACS on the northern side of the 
disk just outside the main dusty ring (Schneider et al. 
2006). This diffuse emission may trace small grains that 
must be replenished because they are being radiatively 
driven from the system. One possible scenario is that 
the northern (J a -band enhancement is the site of recent 
collisional activity; however, small grain structures are 
expected to disappear on 1000 year timescales either be- 
cause small grains are removed by radiation pressure or 
they are smeared out azimuthally by orbital motion. 

If the dust grain density is higher in the northern arm, 
then previous models of the scattering phase function 
may need to be revised. Schneider et al. (2006) fit the 
1.1 fan NICMOS scattered-light image of the HD 181327 
disk using a Henyey-Greenstein-like phase function with 
an asymmetric scattering parameter, g — 0.30 ± 0.03, 
assuming that the grains are uniformly distributed in az- 
imuth. This model implies that the scattered light from 
the northern 120° wedge is 2.1 ± 0.1 times brighter than 
the scattered light from the southern one. If, however, 
the grain density in the northern 120° wedge is 1.4 ± 0.1 
times larger then that in the southern one, then grains 
on the northern side forward scatter (2.1/1.4=) 1.5x the 
amount of light compared to grains on the southern side, 
corresponding to a revised, g = 0.16±0.04. Since the 
physical origin of the thermal brightness enhancement 
on the northern side is currently unknown, 0.12 < g < 
0.33. 

Our PSF-subtracted Q a -band image shows possible 
emission interior to the 86.3 AU ring; however, this struc- 
ture may have been generated by changes in the PSF 
structure during the 5.5 hours over which our observa- 
tions were made. Warm dust components have been ob- 
served interior to the dusty rings around HR 4796A and 
Fomalhaut. Koerner et al. (1998) and Stapelfeldt et al. 
(2004) have argued that these grain populations may be 
generated by collisions among asteroid-like parent bod- 
ies. High resolution thermal-emission imaging at shorter 
wavelengths is needed to confirm the presence of a pos- 
sible warm dust component around HD 181327 at higher 
spatial resolution. 

4. A TWO GRAIN MODEL 

4.1. Thermal-Emission 

Single-temperature black bodies grains and power-law 
size distributions of olivine grains are inadequate for re- 
producing the observed HD 181327 grain properties. (1) 
Black body grains in radiative equilibrium, located at a 
distance of 86 AU around aL, = 3.1 Lq, T e ff = 6450 
K star (Schenider et al. 2006), are expected to possess a 
grain temperature, T gr = 40 K, significantly cooler than 
that inferred from the infrared SED (80 K). (2) Single 
temperature black body fits to the IRS and MIPS SED- 
mode data are unable to reproduce the overall shape 
of the spectrum. The photosphere-subtracted SED for 
HD 181327 (Fig. 2]) possesses a sharp rise in the flux at 
wavelengths >20 /im that peaks at 60-75 /im. Adequate 
fits to IRS continuum produce only three quarters of the 
observed 60-75 /im flux. (3) Power-law size-distributions 
of olivine grains (with a m i n determined by the minimum 



blow-out size and a max = 20 /im) that successfully re- 
produce the gross features of the HD 181327 SED scatter 
0.2 mJy arcsec -2 of the incident stellar light at 1.1 /im 
and 90° scattering, approximately 25% of that observed. 

While the MIPS SED-mode data have very low resolu- 
tion (R ~ 15-25) and do not cover the 40 - 55 fan range, 
the 60-75 fan peak indicates emission from a relatively 
broad solid state feature. Since water ice is expected to 
be abundant and crystalline water ice possesses a well- 
studied emission plateau at 40 - 70 fim with features 
at 43 and 62 fan (Omont et al. 1990), we model the 
far-infrared SED using amorphous olivine grains to cre- 
ate the underlying continuum and water ice (polluted 
with 0.1% amorphous carbon by volume) to reproduce 
the observed emission features. We calculate absorption 
coefficients using laboratory measurements of the optical 
constants for amorphous olivine (MgFeSiO^ Dorschncr 
et al. 1995), water ice near freezing (Warren 1984), and 
amorphous carbon (Preibisch et al. 1993) using Brugge- 
man effective medium and Mie theory (Bohren & Huff- 
man 1983). 

We model the observed SED by summing over the 
thermal-emission produced by each species, requiring 
that both the amorphous olivine and water-ice grains are 
located at 86.3 AU, consistent with the 1.1 fan scattered- 
light image 
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where a is the grain radius, dn/da cx a p is the grain 
size distribution with minimum grain radius, a min , and 
maximum grain radius, a max ; Qv, a bs,. 



and T s are the 
absorption efficiency and temperature of grain species, s, 
and d is the distance from the observer to the star. We 
estimate the minimum grain size assuming that grains 
with /3 — F ra d/ F grav < 0.5 will be gravitationally bound: 
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(Burns et al. 1979; Artymowicz 1988), where (= 3.1 
Lq) and (= 1.4 M©) are the luminosity and mass of 
HD 181327, (Qpr, s (a)) = (J Q pr {a, \)F x dX) {J FxdX)- 1 
is the radiation pressure efficiency averaged over the stel- 
lar spectrum, and p s is the species density. For HD 
181327, we estimate minimum grains sizes, a m i n = 1.0 
and 1.5 fim for amorphous olivine (p s — 3.71 g cm~ 3 ) 
and water ice (p s = 0.92 g cm -3 ), respectively. We set 
amax = 20 /im. Since 2ira max > 100 /im, the presence 
of larger grains should not dramatically impact our SED 
model. We let p be a free parameter and require that p 
be the same for both the olivine and ice distributions. 

We reproduce the observed SED using 6.4xl0 25 g 
amorphous olivine grains with dn/da oc a~ 2,25 (with es- 
timated T gr = 48 - 93 K) and 2.0xl0 25 g of water ice 
grains also with dn/da oc a~ 2 25 (with estimated T gr — 
39 - 50 K; see Fig. [4]). Our model SED possesses a re- 
duced \ 2 value of 22; it is slightly deficient in flux at 
wavelengths 30 fan < A < 55 fan. Our T-ReCS image 
raises the possibility that there may be additional dust 
grains interior to the ring that we have not included in 
our model. If such warm dust grains exist, then their 
emission may account for the discrepancy between our 
observed and synthetic SEDs. Numerical simulations of 
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collisions among dust grains around early-type stars sug- 
gest that radiation pressure modifies the size distribu- 
tion from theoretical collisional equilibrium by removing 
grains with sizes just above the blow-out size (Thebault 
& Augereau 2007; Krivov et al. 2006). The shallow grain 
size distribution inferred from our modeling may indi- 
cate that small grains adjacent to the blow-out size have 
been removed by radiation pressure. Indices of refraction 
have been measured for water ice created at temperatures 
lower than freezing in which the ice is more amorphous 
(Hudgins et al. 1993), and the 42 and 63 /im emission 
features for 100 K amorphous water ice are significantly 
less pronounced than for water ice near freezing. The 
overall better fit provided by crystalline water ice may 
indicate that the ice has been heated in the past and has 
since cooled to the present temperature. 

4.2. Scattered- Light 

We estimate the scattered-light surface brightness, 
generated by our grain population, with scattering ef- 
ficiency, Q,v,sca,s, radius, a, and size column density dis- 
tribution, dNg r ^ s /da, 

SB = Y f ^l fs(0)Q ms na 2 ^^da (3) 

(Augereau & Bcust 2006), where L, )V is the specific lu- 
minosity of the star, D is the distance to the grains, and 
f s (9) is the scattering phase function of species, s. The 
specific luminosity of the star can be inferred from its 
observed flux and distance, >v = kndrF*^, where d = 
50.6pc is the stellar distance and F*^(l.l/xm) = 5.9 Jy 
(Schneider et al. 2006) is the inferred stellar flux in the 
NICMOS filter. The scattering phase function can be 
approximated by the Henyey-Greenstein phase function, 
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(Henyey & Greenstein 1941), where the scattering angle, 
0, is the angle of deviation from forward scattering. If 
the dust is azimuthally symmetric, then we estimate a 
1.1 /Ltm scattered-light surface brightness, SB = 0.81 mJy 
arcsec -2 at D a = 86.3 AU and = 90° (the disk ansa), 
consistent with observations by Schneider et al. (private 
communication) . 

In addition to the scattered-light surface brightness, 
the ACS and NICMOS observations make two addi- 
tional measurements of the grain population around HD 
181327. Schneider et al. (2006) report a scattered-light 
disk color [F606W]-[F110W] = 0.5 ± 0.3 mag at D a = 
86.3AU. The dust grains in our model possess a color 
[F606W]-[F110W] = -0.07 mag, 1.9a lower. Schneider et 
al. (2006) also report an asymmetric scattering factor 
g = 0.30±0.03 which we revise to 0.12 < g HG < 0.33 
in our discussion of the resolved thermal-emission (§3). 
The dust grains in our model possess (g) = 0.87, sub- 
stantially higher than inferred from the scattered light 
observations, because water ice forward-scatters light ef- 
ficiently independent of grain size. The large discrepancy 
between the observed and inferred asymmetric scattering 
parameters may be the result of the unrealistic assump- 
tion of spherical grains. Models including more realistic 
grain shapes may decrease the inferred value of g while 
reproducing the overall thermal emission and scattered 



light from the disk. More detailed modeling of the dust 
grains in this system is needed to reproduce all of the 
observations. 

5. DISCUSSION 

Recent theoretical work suggests that /zm-sized water 
ice grains around A- and F-type main sequence stars are 
more effectively destroyed via photo-desorption by stellar 
ultra-violet photons than by sublimation (Grigorieva et 
al. 2007). We estimate the water ice grain erosion rate 
due to photo-desorption 

a = —r 1 ( 5 

where r\ is the fraction of the surface covered by ice, 
m (H20) is the mass of a water molecule (= 3xl0~ 23 g), 
Pi is the density of water ice (=0.92 g cm -3 ), Y is the 
desorption efficiency at ultraviolet wavelengths (=10 -3 ), 
and N a bs is the flux of absorbed photons at the location 
of the water ice grains, 

Nabs = J Xmin A^kc Qabs{X)dX ' (6) 

where L\ is the specific stellar luminosity, \ m in (= 0.091 
jtxm), and \ max (= 0.24 ^m) (Grigorieva et al. 2007). 
We estimate that the flux of absorbed ultra-violet pho- 
tons that desorb water molecules from ice grains, N aos — 
4.3 xlO 11 photon s _1 cm~ 2 , at D Q = 86.3 AU, assuming 
HD 181327 is 10 x brighter at ultra-violet wavelengths 
than inferred from the stellar photosphere model, con- 
sistent with TD1 observations (Thompson et al. 1978). 
We calculate a = 1.1 xl0~ 3 /im yr -1 if the surface of 
the grains are 99.9% water ice and 0.1% amorphous car- 
bon (77 = 0.999), corresponding to a photo-desorption 
lifetime of 1400 yr for grains with radii, a = 1.5 /j,m, 
similar to the estimated collision timescale for grains in 
this system (^1300 yr; Chen et al. 2006). Since the 
photo-desorption timescale is significantly shorter than 
the stellar age (~12 Myr), we hypothesize that icy grains 
are replenished through collisions among icy parent bod- 
ies, similar to Kuiper Belt objects in our solar system. 

The cold grain temperatures (T gr < 100 K; Meyer et 
al. 2007) in debris disks has led to speculation that these 
systems are icy Kuiper Belt analogs; however, the dif- 
ficulty in detecting mid-infrared spectral features from 
these systems (Chen et al. 2006; Jura et al. 2004) makes 
confirmation of this supposition challenging. To date, the 
most evolved young, circumstellar disk with observed wa- 
ter ice absorption/emission features is the isolated, Her- 
big F7IIIe star HD 142527. ISO SWS and LWS spectra 
of this source detect both the 43 and 62 /im emission 
features associated with the presence of cold, predomi- 
nantly crystalline, water ice (Malfait et al. 1999). De- 
bris disks around main sequence stars are fainter than 
their younger Herbig Fe counterparts, making detection 
of water ice emission even more challenging. The possi- 
ble evidence indicating the presence of water ice emission 
around HD 181327 is the first for a debris disk. Approx- 
imately 25% of the cold dust around HD 181327 may be 
water ice, significantly less than the 60% observed toward 
HD 142527. The relative deficiency of water ice around 
HD 181327 may indicate that it has been incorporated 
into larger bodies (that are less efficient emitters) and/or 
removed via photo-desorption and radiation pressure. 
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6. CONCLUSIONS 

We have obtained a T-RcCS Q a -band (18.3 /um) im- 
age and a 55 - 90 /im MIPS SED-mode spectrum of 
HD 181327, an F5/F6V star in the (3 Pic moving group 
with a resolved scattered-light disk. Based on our anal- 
ysis of the multi-wavelength imaging and spectroscopy, 
we conclude the following: 

1. The northern arm of the HD 181327 disk is 1.4 times 
brighter than the southern arm in thermal-emission, sug- 
gesting that either the density and/or the temperature 
of the grains in this region is higher than in the southern 
arm of the disk. 

2. The overall properties of the unresolved spectral en- 
ergy distribution (at 5.5 - 90 M m )j the resolved thermal- 
emission at 18.3 /im, and the resolved scattered-light at 
1.1 /im can be reproduced by a population of 1 - 20 /im 
amorphous olivine and crystalline water ice grains lo- 
cated at a distance of 86.3 AU from the central star. 

3. Since the estimated photo-desorption lifetime of 1.5 
jivci water ice grains is 1400 yr, significantly shorter than 
the age of HD 181327, the grains must be replenished 
from a reservoir such as collisions among parent bodies, 



perhaps Kuiper Belt objects. 
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Fig. 1. — (a) Q a -band observation of HD 181327. Contours are spaced according to the smoothed 1-<t background noise level, (b) PSF 
and photometric reference star a Ind, with contours scaled to match the brightness levels in panel (a). Both images' contours have been 
smoothed by a FWHM=CK'5 Gaussian. 
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Fig. 2. — (a) Enclosed flux curves for circular and elliptical (e = 0.525, P.A.= 107°) apertures. The curve labeled "total" corresponds to 
the image in Fig.[T^,, while the "reference" curve is derived from the PSF star a Ind, scaled to the estimated photospheric flux (F u> * = 53 mjy) 
of HD 181327. The "subtracted" curve is the difference between the two, showing a detected Q a -band excess of 61.5±2.7mjy inside a 3" 
circular aperture, (b) Azimuthally averaged brightness profiles derived from the enclosed flux curves for HD 181327 and a Ind (scaled to 
Fu,*)- (c) Same as (b) from the PSF-subtracted curve in panel (a). AH panel s show vertical lines marking e- folding distances which bracket 
the peak of the Z/5T/NICMOS scattered-light profile (Schneider ct al. 200(J). Error bars are derived from background noise levels, and do 
not include systematic errors in the flux zero point, and subtracted profiles do not include PSF uncertainty. 
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Fig. 3. — Residual image of HD 181327 after subtraction of the stellar PSF. The image here has been smoothed by a FWHM=0'/5 
Gaussian. Contours are spaced according to the 1-cr smoothed background noise level. The magenta/green tra nsition occurs at the 2-0- 
level. Ellipses show the e-folding distances which bracket the peak of the ffST/NICMOS scattered-light profile (Schneider ct al. 2006). 
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Fig. 4. — Photosphere-subtracted IRS and SED-mode spectra for HD181327, shown with black error bars. Overlaid are T-ReCS (this 
work), MIPS (Smith et al. 2006), and IRAS photometry (not color-corrected), shown with red, green, and blue error bars. Model emission 
from amorphous olivine and crystalline water ice are shown with dashed and dotted lines, respectively. The final model for the system is 
shown with a solid black line. 



